from Excitation_Carbon import *

cs = CS_Excitation(species_name = "C+4")


id = 0

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =299.0,
    A   = -1.820e-5,
    B   = 1.093e-3,
    C   = 1.675e-2,
    D   = -8.788e-3,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+4_1s21S-1s2s3S",
    X1  = 1.02,
    P   = -1.777,
    Q   = 1.826)


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =304.0,
    A   = 2.830e-2,
    B   = -3.159e-2,
    C   = 2.694e-2,
    D   = -1.063e-2,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+4_1s21S-1s2s1S")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 304.0,
    A   = 8.926e-5,
    B   = -2.824e-3,
    C   = 7.379e-2,
    D   = -2.670e-2,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+4_1s21S-1s2p3P")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =308.0,
    A   = 1.642e-2,
    B   = -7.345e-2,
    C   = 9.320e-2,
    D   = 0.0,
    E   = 1.247e-1,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+4_1s21S-1s2p1P")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =354.0,
    A   = 5.061e-5,
    B   = -5.619e-4,
    C   = 2.116e-3,
    D   = 7.710e-4,
    E   = 0.0,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+4_1s21S-1s3d3D")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =355.0,
    A   = 1.286e-2,
    B   = -2.133e-2,
    C   = 1.656e-2,
    D   = 0.0,
    E   = 2.107e-2,
    F   = 0.0,
    omega = 1.0,
    equation_number = 6,
    reaction_name = "Excitation_C+4_1s21S-1s3p1P")


cs.plot_cs()
cs.save_cs()
cs.cal_total_cs()